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Abstract 

A family of spherical shells with varying thickness is derived by 
using a simple Newtonian potential-density pair. Then, a particular 
isotropic form of a metric in spherical coordinates is used to construct 
a General Relativistic version of the Newtonian family of shells. The 
matter of these relativistic shells presents equal azimuthal and polar 
pressures, while the radial pressure is a constant times the tangential 
pressure. We also make a first study of stability of both the Newtonian 
and relativistic families of shells. 

Keywords: gravitation. 

1 Introduction 

Distributions of matter in form of shells have been a useful tool in astro- 
physics and General Relativity, often providing simplified but analytically 
tractable models in cosmology, gravitational collapse and supernovae. In 
General Relativity, the pioneering work on constructing solutions of the 
Einstein field equations representing shells was done by Israel [I]. In his 
formalism, the shells are constructed by joining two manifolds across a 
three-dimensional time-like surface, and the matter representing the shell 
is interpreted as a discontinuity of the extrinsic curvature on the surface. 
This formalism has been used to study several models of relativistic shells. 
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[2] give a brief review of some of these models and also discuss shells which 
admit closed equations of state, including non-linear polytropic models. 

An exact solution of a static shell surrounding a black hole was obtained 
by [3] , and the stability of such a configuration against radial perturbations 
was studied by [3], while non-radial perturbations were analysed by [5]. Non- 
radial perturbations of self-gravitating static fluid shells in the Newtonian 
gravity were studied by [6]. 

In this work, we propose a Newtonian analytical potential-density pair 
that represents a family of spherical shells with varying thickness. For this 
potential, we calculate the rotation curves of test particles and make a first 
study of stability of the shells by considering the stability of circular orbits 
of the test particles. This is done in Section [2j In Section [3l we consider 
a General Relativistic version of the Newtonian shells. For this task, we 
employ a particular isotropic form of a metric in spherical coordinates. We 
calculate the components of the energy-momentum tensor of the shell's mat- 
ter and also analyse the geodesic circular motion of test particles and discuss 
their stability. Finally, in Section 0] we summarize our results. 



2 A Model of a Newtonian Shell 

In this section, we present a simple potential-density pair that describes a 
family of spherical shells of matter. We shall consider the following potential: 

(r n + b n ) 1 

where b > is a parameter with dimensions of length and n is a non-negative 
integer. By using the Poisson equation in spherical coordinates, 

1 d / 9 d$\ 

V— =4vrG/9, (2) 



(3) 



r 2 dr V, dr J 
one obtains a mass density given by 

m(n + l)b n r n ~ 2 



Air (r n + b n ) 



,7^2+l/rT 



For n = 2, the pair (fTJ) and ([3]) is reduced to the Plummer model [3(8] that 
has a monotone decreasing density profile. However, for n > 2 the density 
vanishes on r = and thus we have a shell-like distribution of matter. Fig. 
[TJa) shows some curves of the dimensionless mass density p = p/(m/b 3 ) as 
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Figure 1: (a) The mass density p = p/(m/b 3 ) (equation [3]) , as function of 
r/b for n = 3, n = 6 and n = 9. (b) The circular velocity v c = v c / (Gm/b) 1 ! 2 
(equation [5]) , and (c) the angular momentum h = h/(Gmb) 1 ^ 2 (equation [6]) 
as functions of r/b for n = 3, n = 6 and n = 9. 



function of r/b for n = 3, n = 6 and n = 9. For higher values of n, the mass 
density becomes more concentrated (thinner shell), with maximum at 



n 



n + 3 



l/n 



b. 



(4) 



For large values of r, the density (JTJ) decays as l/r n+3 , thus very fast for 
large n, as can be seen in Fig. [T^a). 

The circular velocity v c = (r^ ,,) 1 / 2 and angular momentum per unit 
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mass h = rv c of test particles moving in circular orbits are given by 



/ Gmr n 

" C = V(^ + 6«) 1+1/ "' (5) 
/ Gmr n + 2 

V(r n + & n ) 1+1/n ' 

Curves of the circular velocity v c = v c /(Gm/b) 1 ^ 2 and of the angular mo- 
mentum h = h/(Gmb) 1 ^ 2 for n = 3, n = 6 and n = 9 are depicted in Figs 
HKb) and (c), respectively. The stability of circular orbits can be determined 
using the Rayleigh criterion of stability for a rotating fluid ([9]; see also 
[TU]). and that can be adapted to study circular orbits of test particles in 
self-gravitating systems (see e. g. for an application of the Rayleigh cri- 
terion of stability in General Relativity). The condition for a stable orbit 
is 

4 > 0. (7) 
dr 

Using equation ([6j) , it is straigthforward to check that the orbits are always 
stable. 

It is instructive to consider the limit n — > oo. In this case, we have an 
infinitesimal thin shell on r = b, and the limits of the potential ([T|) and 
circular velocity jSJ are 

$ = { n ~ (8) 

if r > &> w 

[0 if < r < b, 

VC= \^ ifr>6- (9) 

Thus, for r < b the potential is constant and a test particle is at rest inside 
the shell, and if r > b the shell exerts on the test particle the same gravita- 
tional force as a point mass located on r = 0, in agreement with Newton's 
Shell Theorem. 



3 A Model of a General Relativistic Shell 

Now, we study a General Relativistic version of the Newtonian potential- 
density pair discussed in the previous section. For this, it is convenient 
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to choose a particular metric in isotropic form in spherical coordinates 
(t, r, 9, (p), which we write as 



ds 2 = (j-^jj 2 c 2 dt 2 - (1 + /) 4 (dr 2 + r 2 d9 2 + r 2 sin 2 6dp 2 ) , 



(10) 



where / = f(r). In this particular form of metric, the Schwarzschild solution 
is given by / = Gm/(2c 2 r). A similar form of the metric (llOj) in cylindrical 
coordinates was used by [12j to construct General Relativistic versions of 
Newtonian models for distributions of matter in galaxies. 

For metric (fTU|) . the Einstein equations = — (87rG/c 4 )T^2, yield the 
following expressions for the non-zero components of the energy-momentum 



tensor T, 



1 d 



>d/ 

2vrG(l + /) 5 r 2 dr V' dr 



2vrG(l + /) 5 (l-/) dr 

„4 



d/ ff_ + d/ 
r dr 



TS 



rpip 
<P 



4ttG(1 + /)°(1-/) 



/ d 2 / + /d/ 
dr 2 r dr 



dr 



(11) 
(12) 

(13) 



Because T^ v has only diagonal non-zero components, the energy density is 
directly given by e = T*/c 2 and the pressures or tensions along a direction k 
read P k = 
P<p/c 2 reads 



T%. The 'effective Newtonian density' pn = e + P r /c 2 + P s /c 2 + 



1 d 



Pn 



2ttG(1 + /) 5 (1-/) r 2 dr 



-d/ 
dr 



(14) 



The relation between the function /(r) and the Newtonian potential $(r) 
is obtained by comparing equation (fT4l) or (fTT|) with equation ([2]) in the 
non-relativistic limit when 



2?' 



/ ■- — ■ (15) 
Now, using equation (fT5|) with the potential ([I]), equations (fTT]) - (fT4l) can 
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be cast as 



(n + l)b n r 



(16) 



e = 



5 - 



4^2-4/n (1 



n6 n f n ~ 2 



(17) 



-Pfl — -Pos — 



8vr £2-4/« (_! + £l/n) ^ + £l/„) 
(n + l)6 n f n " 2 



5 ' 



(18) 



PN = 



4^2-5/n (_! + £l/n) ( X + ' 



(19) 



where the dimensionless variables and parameters are as follows: f = r/r s , 
b = b/r s , e = e/ (m/rf ), pjv = Pn / {m/rf)-, Pk = Pk/(mc 2 As)> and where we 
defined r s = Gm/(2c 2 ) and £ = f n + 6 n - We assume again that n > 2. For 
n = 2, equations (fT6]) - (fT9|) reduce to the perfect fluid sphere derived by [13], 
which has as equation of state similar to the classical polytrope of index 5. 
Note that the relation between P r and Pg is independent of the radius and 
is simply P r = 2Pg/n. 

In order to be physically meaningful, the components of the energy- 
momentum tensor should satisfy the energy conditions. The strong energy 
condition states that pn > 0, whereas the weak energy condition imposes 
the condition e > 0. The dominant energy condition requires |-Pfc/ e l < c 2 . 
Equations fjl6(> — (119f) show that the weak energy condition is always satisfied, 
and the strong condition holds if b > 1. This condition also ensures that 
we have pressures everywhere. The inequality |-P r /e| < c 2 is satisfied when 
b > (n + 2)/(n + l), and we have \P e /e\ < c 2 if b> (3n + 2)/(2n + 2). Thus, 
there exists an interval of the parameter b where all energy conditions are 



Figs [2ja)-(d) display, respectively, some curves of the energy density 
e, the effective Newtonian density pjy, the radial pressure P r and the az- 
imuthal/polar pressure Pg as functions of the radial coordinate for the shell 
with n = 3 and parameter values 6 = 2, 2.5 and 3. In Figs EJa)-(d), the 
same quantities are displayed for the shell with n = 6 and parameter values 
b = 2.5, 3 and 3.5. All quantities have the profile of a shell-like distribu- 
tion of matter. As the value of the parameter b is increased, the densities 
and pressures are smoothed out. The densities and pressures become more 
concentrated when n is increased. 

As we did in the Newtonian model, we also study the circular orbits and 
stability of test particles. Without loss of generatity, if the geodesic motion 



satisfied. 
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Figure 2: (a) The energy density e = e/(m/rf) (equation [T6|) . (b) the ef- 
fective Newtonian density pN = pn / {fn/r^) (equation [T9j) . (c) the radial 
pressure P r = P r /(mc 2 /H?) (equation 1 1 7j) and (d) the azimuthal/polar pres- 
sure Pq = Pg/(mc 2 /rg) (equation 1 18|) as functions of f = r/r s for n = 3, 
and parameter b = b/r s = 2, 2.5 and 3. 



7 



2 
1.5 
e 1 
0.5 





x 10" 



(a) 





— b=2.5 

- - b=3.0 

b=3.5 


/ ' A* v 
/ ' ' \ \* 
/ ' \ \ 









x 10" 



1.5 



0.5 



(c) 



2 4 
r 





— b=2.5 

- - b=3.0 

b=3.5 


/ / V 

/ / • ' \ ^ * 
/ 1 ■' \ v \ " 






Figure 3: (a) The energy density e = e/(m/rf) (equation [T6|) . (b) the ef- 
fective Newtonian density pN = pn / {fn/r^) (equation [T9j) . (c) the radial 
pressure P r = P r /(mc 2 /H?) (equation 1 1 7j) and (d) the azimuthal/polar pres- 
sure Pq = Pg/(mc 2 /rg) (equation 1 18|) as functions of f = r/r s for n = 6, 
and parameter b = b/r s = 2.5, 3 and 3.5. 
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n 3456789 
b^n. 1-78 2.00 2.18 2.32 2.43 2.52 2.60 

Table 1: Approximate minimum values of the parameter b = b/r s that 
ensures sublumial rotation curves for some values of n. 

is confined on the = tt/2 plane, the expressions for the circular veloc- 
ity and for the angular momentum are formally the same as in cylindrical 
coordinates, and we have P21 



(20) 



»-* '^fr ^ir-p + 2 J 4(w)r (21) 




For the shell potential J!]), we get 



y (-1 + £}l n ) (b n - f n + £ 1+ V n ) ' 

^3/(2n) ^1 fn \jn _|_ £ 1+2/n ^n^l/n ' 



(22) 



(23) 



where v c = v c /c, h = r s c and £ = f n + b n . 

Some curves of the circular velocity v c = v c /c and angular momentum 
h = r s c for relativistic shells are shown in Figs IU(a)-(d). In Figs H^a)-(b), 
the curves of rotation are displayed for the shells with n = 3 and n = 6, 
respectively, and the corresponding curves of angular momentum are plotted 
in Figs H|c)-(d). As the values of the parameter b = b/r s is increased, the 
velocities becomes less relativistic. The minimum value of b such that the 
rotation curve is sublumial has to be found numerically. Table [T] shows some 
approximate minimum values for b for shells with n = 3 to n = 9. Note that 
these minimum values are always larger than those required by the energy 
conditions. 

The curves of angular momentum in Figs[4]Jc)-(d) show regions of insta- 
bility at high relativistic velocities. In the Newtonian limit, the orbits are 
always stable, thus this instability is a purely relativistic effect. Using the 
Rayleigh criterion of stability, curves of dh/df = as functions of f and b 
are shown in Fig. [5] for models of shells with n from n = 3 to n = 6. In each 
graph, the minimum value of b is that given in Table HJ In the inner part of 
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Figure 4: (a)-(b) The circular velocity v c = v c /c (equation |22|) . as a function 
of r = r/r s for the relativistic shell model with (a) n = 3 and (b) n = 6. 
(c)-(d) The angular momentum h = r s c (equation [23]) . as a function of 
f = r/r s for the relativistic shell model with (c) n = 3 and (d) n = 6. 
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the curves, orbits are unstable. The regions of instability become larger as 
the parameter n is increased. 

4 Discussion 

We proposed a very simple potential-density pair in Newtonian gravity that 
describes a family of spherical shells whose thickness is function of a param- 
eter n. In the limit n — > oo, one obtains an infinitesimal thin shell. This 
family of shells is found to be stable by the Rayleigh criterion of stabil- 
ity. We also studied a General Relativistic version of the Newtonian model 
of shells by using a particular metric in isotropic form in spherical coordi- 
nates. These relativistic shells have radial pressures that are different from 
the equal azimuthal and polar pressures, but are proportional to for a fixed 
value of n. Also, parameters can be chosen such that the matter in the 
shells satisfies all the energy conditions. The Rayleigh criterion of stability 
shows that the relativistic shells have regions of unstable circular orbits of 
test particles, and this instability is a pure relativistic effect. 
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